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This article describes the application of recently introduced complex networks concepts and meth- 
ods to the characterization and analysis of cortical bone structure. Three-dimensional reconstruc- 
tions of the system of channels underlying bone structure are obtained by using histological and 
computer graphics methods and then represented in terms of complex networks. Confluences of 
two or more channels are represented as nodes, while the interconnecting channels are expressed as 
edges. The hierarchical backbone (the tree with maximum depth) of such a network is obtained and 
understood to correspond to the main structure underlying the channel system. The remainder of 
the network is shown to correspond to geographical communities, suggesting that the bone channel 
structure involves a number of regular communities appended along the hierarchical backbone. It 
is shown that such additional edges play a crucial role in enhancing the network resilience and in 
reducing the shortest paths in both topology and geometry. The recently introduced concept of frac- 
tal dimension of a network ( cond- mat/0503078 ) is then correlated with the resilience of the several 
components of the bone channel structure to obstruction and failure, with important implications 
for the understanding of the organization and robustness of cortical bone structure. 



I. INTRODUCTION 

Most biological structures and phenomena involve in- 
teractions between several components, required in or- 
der to obtain proper functionality and behavior. For 
instance, the nervous system involves a large number 
of neuronal cells interconnected through synapses, while 
the skeleton consists of many bones attached one an- 
other. Because of such a type of organization, biological 
structures are often properly represented and modeled 
in terms of graphs or networks of interconnected nodes. 
Therefore, the nervous system can be modeled in terms 
of a graph whose nodes are understood to represent the 
neurons and the edges the synaptic connections (e.g. 0). 
Interestingly, despite the continuing interest in complex 
networks research (e.g. 0,0), relatively little attention 
has been focused on their use to represent geometrical 
structures in biology. One of the few exceptions is the ap- 
plication of complex networks to analyze the channel net- 
work underlying cortical bones, which was reported only 
recently Q . The current work considerably extends that 
preliminary work in several important aspects by taking 
into account recently introduced powerful concepts |5j 
such as the fractal dimension of a network |fj , the hierar- 
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chical backbones of networks |jj , as well as communities 
(e.g. @), resilience to attack (e.g. 0) and shortest paths. 
By treating such concepts in an integrated fashion, it has 
been possible to infer a series of insights about the net- 
work of channels underlying cortical bone structures. 

II. THE HAVERSIAN SYSTEM 

Visual observation of the diaphysis of transversally sec- 
tioned long bones immediately reveals that this part in 
the bone is formed mainly by compact structures. How- 
ever, at microscopic level, such a structure is verified to 
be composed of collagen fibers organized as lamellas with 
3 to 7 /im of thickness which are either parallel one an- 
other or concentric to the microscopy canal, constituting 
the Havers System or Osteon. Such a network, which 
is called the Havers channels, follows the long axis of 
the bone, communicating with itself, with the medullary 
cavity and with the external surface of the bone through 
transversal or oblique channels called in honor of Volk- 
mann. These channels have as main function to contain 
the sanguineous vases to nourish and supply oxygen to 
the deepest cells of the bone tissue. During bone growth, 
its structure is constantly modified by undoing old and 
forming new Havers systems. In case the distribution of 
blood vessels through the bone were implemented exclu- 
sively by a dichotomic tree, the elimination of specific 
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FIG. 1: Non geographical Haversian Network (a) and the 
longest extracted tree, understood as the respective hierar- 
chical backbone (b). 



channels could hinder nutrient and oxygen supply. Some 
type of architectural redundancy is therefore expected 
in order to compensate for eventual channel elimination 
through anastomosis during the growth or pathologies in 
adults. 



III. METHODOLOGY 

The umerous of the adult cat used in this work was 
obtained from an animal necropsed in the Pathology 
Sector of the Veterinary Hospital of the Federal Uni- 
versity at Uberlandia, Brazil. After dissection, a 5cm 
thick bone ring was extracted and sliced according to 
histology traditional procedures. After digitally imaged 
(700 x 800 pixes), the sections were registered in size, 
translation and orientation, and the Havers and Volk- 
man channels were manually identified, in order to allow 
three-dimensional reconstruction of the Haversian system 
by using computer graphics techniques |lO| . A undirected 
graph was extracted from the three-dimensional struc- 
ture by assigning a node to each channel confluence and 
an edge to each channel. Channels shaped as a "V" were 
understood to correspond to confluence of two channels, 
being therefore represented as a node with degree 2. The 
complete graph included 852 vertices and 1016 edges, im- 
plying average degree (fc) = 2.4. 

Figure ^a) shows the non-geographic representation 
of the network obtained from the Haversian system. The 
respective hierarchical backbone corresponding to the 
longest underlying tree was obtained as described in 
and is shown in Figure ^b). More specifically, a tree was 
extracted starting at each of the network nodes. Such 
trees were obtained by following the hierarchical levels 
from the starting node, while removing all cycles. The 
hierarchical backbone was chosen as the longest of such 
trees, but the widest tree was also identified and consid- 
ered in our investigations. Such an approach was mo- 
tivated by the assumption that the Havers channels are 
meant for blood supply along the bone structure. 

Given the spatial adjacency restrictions imposed by the 
3D space where the network lives pl|. it is reasonable 



to expect a branched channel organization such as those 
found in body vascularization and lungs. However, by be- 
ing internal to a rigid structure (bones), such a structure 
would imply a serious shortcoming in the sense that even- 
tual obstructions of one of the channels would mean the 
congestion of all downstream channels Therefore, 
we expect to find additional channels providing bypasses 
along the bone channel system, yielding a more resilient 
hybrid network. 

In order to verify the hypothesis above, the original 
channel network was removed from the original Haver- 
sian network (by using the concept of difference between 
two graphs). The obtained results are shown in Figure|21 
which includes the two-dimensional projection of the ge- 
ographic bone network (center) with its five identified ge- 
ographical communities (different gray levels), as well as 
the respectively subnetworks defined by each geograph- 
ical community. Note that the edges shown in the cen- 
tral structure in the previous figure correspond to those 
which are part of cycles in the original net. Such a fact 
is compatible with the above hypotheses that the Haver- 
sian system involves a tree backbone complemented by 
additional bypass edges (those in the five communities), 
also indicating that the latter are organized so as to form 
well-defined geographical communities. Interestingly, the 
same communities were obtained when considering the 
widest hierarchical backbone. 

The loglog distributions of the node degrees of the 
Haversian network and its tree backbone are shown in 
Figure E^a), from which it is clear that both graphs ex- 
hibit quite similar node degrees. The node degree distri- 
butions of regular, random and scale free equivalent net- 
works, i.e. with the same number of edges and nodes, are 
shown in Figure EJb-d), clearly illustrating that the node 
degree structure of the Haversian system is not similar 
to any of these alternative models, except the scale-free 
case in (d). 

Because the resilience of the Haversian network is one 
of the most important features from the functional point 
of view, we subjected the original bone network to ran- 
dom edge attack, and the result (size of dominant cluster) 
is shown in Figure^] (triangles) in terms of the fraction of 
removed edges. This figure also includes results obtained 
for equivalent tree, regular, random and scale- free models 
with approximately the same number of nodes and edges. 
Note that the Haversian network presents resilience be- 
havior which is quite similar to that of regular networks, 
suggesting that the bone structure could have a marked 
degree of regularity. However, the node degree distribu- 
tions shown in the loglog diagrams in Figure [3] clearly 
indicate that this is not the case. In other words, the 
investigated Haversian channel system presents a node 
degree distribution which is similar to its tree backbone 
(see Figure |3Ja)) or even a scale- free model, but a sub- 
stantially enhanced resilience which follows closely that 
of the equivalent regular model. The explanation for such 



FIG. 2: Haversian Network partitioned into five well-defined geographical communities (center) and the respectively defined 
subgraphs. 
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gree of topological self-similarity of a network. Several 
interesting networks, including protein interaction and 
internet, have been verified to exhibit such a fractal or- 
ganization. The algorithm involves the coverage of the 
network by boxes of size I b , such that all vertices inside 
one of such boxes connect one another at a distance which 
is smaller or equal to Ij, — 1. The number of boxes re- 
quired to cover the whole network, expressed as Nb(Ib), 
is calculated for different values of Is- The topological 
self-similarity is characterized by the following power law 
relationship 
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FIG. 3: Degree distributions of the Haversian system and 
the respective tree backbone (a). The degree distribution of 
regular, random and scale free models with the same number 
of nodes and approximately the same number of edges are 
shown in (b-d). 



mixed topological features is provided by the fact that 
the Haversian system can be understood as a hybrid, in- 
volving a tree backbone to which spatial communities are 
attached in order to provide robustness to flow obstruc- 
tion, as illustrated in Figure |2] 

In order to obtain a more comprehensive characteriza- 
tion of the topological features of the Haversian network, 
we calculated the network fractal dimension as proposed 
by Song, Havlin and Makse |6|, which expresses the de- 



Figures OJa-e) show the loglog curve of Nb/Nt, where 
Nt is the number of nodes in each case, obtained for each 
of the subgraphs, while Figure (5ff) shows the fractal di- 
mension for the complete Haversian network as well as 
its respective tree backbone and the geometrical fractal 
dimension of the bone structure, calculated by the box- 
counting methodology (e.g. ^3)- The topological fractal 
dimensions of equivalent regular, random and scale free 
models are presented in Figure HO The fractal nature 
of the regular network is clear from (a) , while the high- 
est topological fractal dimension is obtained for the scale 
free case in (c). The topological fractal dimension of the 
instances of the attacked Haversian system shown by tri- 
angles in Figure [5] are presented in terms of the rate of 
edge attack in (d), which indicates that the topological 
fractal dimension is negatively correlated with the edge 
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FIG. 4: Despercolation of different models of complex net- 
works 



FIG. 6: Fractal dimension of: (a) regular, (b) random and 
(c) scale-free network. In (d) shows with fractal dimension of 
Haversian network vary in despercolation process. 
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FIG. 5: Fractal dimension of five communities (a-e) and of 
the main Haversian network (f). The fractal dimension of 
the respective tree backbone and the box-counting fractal di- 
mension of the geometrical channel system are also shown in 

(f) 



attack rate, suggesting the use of the former as an indi- 
cation of the resilience of the investigated network. The 
shortest paths between any two nodes in the Haversian 
network and its related tree backbone have also been es- 
timated and the respective histograms are shown in Fig- 
ure[7fa), while the histogram in (b) indicates the shortest 
Euclidean distances. It is clear from these results that the 
addition of the communities to the tree backbone have 
the major effect of reducing substantially the shortest 
paths. 



FIG. 7: Histograms of shortest topological (a) and Euclidian 
distance (b) paths for the Haversian network and its respec- 
tive tree backbone. 



IV. CONCLUDING REMARKS 

The present work investigated the possibility of using 
recently introduced complex networks concepts to rep- 
resent and characterize biological 3D structures, namely 
the Haversian channel system underlying cortical bones. 
The Haversian system of channels has been found to in- 
volve a hierarchical backbone (the longest tree) to which 
a series of geographically well defined communities are 
attached in order to provide bypasses ensuring some de- 
gree of robustness to eventual channel obstruction. Such 
a feature of the Haversian network has been corroborated 
by simulated random edge attack, which indicated that, 
though exhibiting a node degree distribution similar to 
that of an equivalent tree, the bone structure presents 
robustness close to that of a regular network of similar 
size and average degree. Additional quantification of the 
topological redundancy of the Haversian system was ob- 
tained by using the recently introduced concept of net- 
work fractal dimension, which was found to be negatively 
correlated to the intensity of the attack, suggesting the 
use of this measurement as an indication of the network 
resilience. The presence of the geographical communi- 
ties along the tree backbone was also found to contribute 
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significantly to reducing the overall shortest paths in the 
channel system both in topological and distance terms. 
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